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We analyze strong field atomic dynamics semiclassically, based on a full time-dependent description 
with the Hermann-Kluk propagator. From the properties of the exact classical trajectories, in 
particular the accumulation of action in time, the prominent features of above threshold ionization 
(ATI) and higher harmonic generation (HHG) are proven to be interference phenomena. They are 
reproduced quantitatively in the semiclassical approximation. Moreover, the behavior of the action 
of the classical trajectories supports the so called strong field approximation which has been devised 
and postulated for strong field dynamics. 
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I. INTRODUCTION 

In the last two decades multiphoton processes have 
been studied intensively, experimentally as well as the- 
oretically. The inherent time-dependent nature of an 
atomic or molecular excitation process induced by a short 
laser pulse renders a theoretical description problem- 
atic in two respects. Firstly, a full quantum calculation 
in three dimensions requires a large computational ef- 
fort. For this reason, quantum calculations have been 
restricted to one active electron in most cases jj],D . Sec- 
ondly, an intuitive understanding for an explicitly time 
dependent process seems to be notoriously difficult, ex- 
emplified by pertinent discussions about stabilization in 
intense laser fields Many studies have been car- 

ried out to gain an intuitive understanding of the two 
most prominent strong field phenomena, namely High 
Harmonic Generation (HHG) and Above Threshold Ion- 
ization (ATI) . In the well established early analytical for- 
mulation by Keldysh, Faisal and Reiss the atomic poten- 
tial is treated as a perturbation for the motion of the 
electron in a strong laser field ||. 

This picture is still used in more recent models, where 
the classical dynamics of the electron in the laser field is 
explicitly considered, e.g. in Corkum's rescattering model 
which can explain the cutoff observed in HHG for linearly 
polarized laser light in one spatial dimension H]. The 
corresponding Hamiltonian reads J§| 

H = H + E f(t) x sin^i + 6) , (1) 

where Ho = ^p 2 + V(x) is the atomic Hamiltonian, f(t) 
is the time-profile of the laser pulse with maximum am- 
plitude E , and u>q is the laser frequency. The interaction 
of the electron with the atom is specified by the potential 
V. 

Lewenstein et al. extended Corkum's rescattering idea 
to a quasiclassical model which contains one (relevant) 



bound state not influenced by the laser field on the one 
hand side and electrons which only feel the laser field on 
the other side This simple model explains qualita- 
tively the features of HHG well. The same is also true 
for an alternative model, where the electron is bound 
by a zero-range potential (nj However, the basic ques- 
tion if and to which extent these multiphoton processes 
can be understood semiclassically, i.e., by interference of 
classical trajectories alone, remains unanswered. It is as- 
tonishing that no direct semiclassical investigation of the 
Hamiltonian Eq. ([!]) has been performed while a number 
of classical as well as quantum calculations for Eq. (|l|) 
have been published. However, only recently, a semiclas- 
sical propagation method has been formulated which can 
be implemented with reasonable numerical effort. This is 
very important for the seemingly simple Hamiltonian Eq. 
(|l|) whose classical dynamics is mixed and in some phase 
space regions highly chaotic which requires efficient com- 
putation to achieve convergence. Equipped with these 
semiclassical tools we have studied multiphoton phenom- 
ena semiclassically in the frame work of Eq. (]l|). In com- 
parison to the exact quantum solution we will work out 
those features of the intense field dynamics that can be 
understood in terms of interference of classical trajecto- 
ries. 

The plan of the paper is as follows. In section II we 
provide the tools for the calculation of a semiclassical, 
time-dependent wavefunction. In section III we discuss 
Above- Threshold-Ionization (ATI) and work out the clas- 
sical quantities which structure semiclassically the rele- 
vant observables. In section IV we use this knowledge for 
the description of Higher-Harmonic-Generation (HHG). 
Section V concludes the paper with a comparison of HHG 
and ATI from a semiclassical perspective and a short 
summary. 
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II. CALCULATION OF THE SEMICLASSICAL 
WAVE FUNCTION 

A (multi-dimensional) wave function ^(x, t) can be 
expressed as 



*(x,i) = f dx'if(x,x',t)*(x'). 
Jo 



(2) 



Here, Vt'(x') is the initial wave function at t = and 
2f(x, x',i) denotes the propagator. We will not use the 
well-known semiclassical Van Vleck-Gutzwiller (VVG) 
propagator which is inconvenient for several reasons. 
Firstly, one has to deal with caustics, i.e. singularities 
of the propagator, and secondly, it is originally formu- 
lated as a boundary value problem. For numerical ap- 
plications much better suited (and for analytical consid- 
erations not worse) is the so called Herman-Kluk (HK) 
propagator which is a uniformized propagator in initial 
value representation [[ll|Jl2| , formulated in phase space, 



iT-(x,x',i) = ^-L_ JJdpdqC flp (t)e^/ h 



flf 7 (x;q(i),p(i)) 5*(x';q,p) (3) 



with 



(x;q,p) = (~) ' exp f-l (x- q) 2 + rP(x- q, 



9j( 
and 

CqpW 
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(4) 



(5) 



Each phase space point (q, p) in the integrand of Eq. 
(||) is the starting point of a classical trajectory with 
action 5 qp (t) . The terms X y in the weight factor C qp (t) 
are the four elements of the monodromy matrix, X y = 
dx t /dy. The square root in Eq. (|5|) has to be calculated 
in such a manner that Cqp(t) is a continuous function of 
t. The integrand in Eq. (M) is - depending on the system 
- highly oscillatory. Although we restrict ourselves to one 
spatial dimension (see Eq. (ffl)) the number of trajectories 
necessary for numerical convergence can reach 10 7 . We 
note in passing that an integration by stationary phase 
approximation over momentum and coordinate variables 
reduces the HK-propagator to the VVG-propagator jl3) . 

In all calculations presented here we have used a Gaus- 
sian wave packet as initial wave function, 
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With this choice, the overlap 



:(x';q,p) V (x')dx' 



(6) 



(7) 



can be calculated analytically and Eq. (g) reads together 
with Eq. (0) 



dpdq e lS «v {t),h 
Cg P (t) g 1 (x; q(t) ,p(t)) / 7/3 (q, p) (8) 



with a = 7 + /3. For all results presented here we have 
taken 7 = (3. 

For comparison with our semiclassical calculations we 
determined the quantum mechanical wave function using 
standard Fast Fourier Transform split operator methods 

El. 



III. ABOVE THRESHOLD IONIZATION 

We start from Eq. ([!]) with 5 = and use a rectangular 
pulse shape f(t) which lasts for 4.25 optical cycles. This 
setting is very similar to the one used in Jig] . 

The energy spectrum of the electrons can be expressed 
by the Fourier transform of the autocorrelation function 
after the pulse, i.e. for times t>tf, 



00 

a{to) = Re J e wt {^{t)\^ f ) dt , 



(9) 



where "f/ = ^itf) is the wave function after the pulse 
and correspondingly 



(10) 



is calculated by propagating fyf for some time with the 
atomic Hamiltonian Ho only after the laser has been 
switched off. 



A. Quantum mechanical and semiclassical spectra 
for ATI 

We will present results for two types of potentials to 
elucidate the dependence of the semiclassical approxima- 
tion on the form of the potential. 



1. Softcore potential 
First we apply the widely used softcore potential 

MM 



V(x) 



1 



\J x 1 + a 



(11) 



with a = 1 and with an ionization potential I p = 0.670 
a.u.. We have checked that the correlation function 
differs little if calculated with the exact ground state 
or with the ground state wave function approximated 
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by the Gaussian of Eq. (||) where (3 = 0.431 a.u. and 
qp = 0. However, the semiclassical calculation is consid- 
erably simplified with a Gaussian as initial state as can 
be seen from Eqs. (||||). Therefore we use this initial 
state and obtain the propagated semiclassical wavefunc- 
tion in the closed form Eq. (||). In Fig. [j]the quantum and 
semiclassical results at a frequency loq — 0.148 a.u. and a 
field strength Eq — 0.15 a.u. are compared. The Keldysh 
parameter has the value 1.14. The quantum mechanical 
calculation (dotted line) shows a typical ATI spectrum. 
Intensity maxima with a separation in energy of Hloq are 
clearly visible. The first maximum has the highest inten- 
sity while the second maximum is suppressed. 
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FIG. 1. Quantum mechanical (dotted line) and semiclas- 
sical (solid line) ATI spectra for the Hamiltonian of Eq. (^) 
with E — 0.15 a.u., luq = 0.148 a.u. and the softcore potential 
Eq. (0. 

The semiclassical result (solid line) is ambiguous: On 
the one hand there are clear ATI maxima with a sepa- 
ration of Tiljq. All peaks but the first one have roughly 
the correct magnitude. Again the second maximum is 
missing. On the other hand we see a constant shift 
(about 0.02a.u.) of the spectrum towards higher energies. 
Therefore, a quantitative semiclassical description is im- 
possible, at least with the present parameters and the 
softcore potential. Next, we will clarify whether the shift 
in the spectrum is an inherent problem of a semiclassical 
ATI calculation or if it can be attributed to properties of 
the softcore potential. 

2. Gaussian potential 

To this end we take a potential which has been used 
to model the "single bound state" situation mentioned in 
the introduction |17j . It is of Gaussian form 

V(x) = -V exp (-ax 2 ) . (12) 

With our choice of parameters Vb = 0.6 a.u. and 
(j = 0.025 a.u., the potential contains six bound states 
and can be approximated, at least in the lower energy 



part, by a harmonic potential for which semiclassical cal- 
culations are exact. Hence, the semiclassical ATI spec- 
trum with this potential should be more accurate if the 
discrepancies in Fig. ^ are due to the potential and not 
due to the laser interaction. The ground state wave func- 
tion itself is again well approximated by the Gaussian 
Eq. (||) with (3 = 0.154 a.u. and qp = 0. The laser has 
a frequency loq = 0.09 a.u., a field strength Eq = 0.049 
a.u., and a pulse duration of 4.25 cycles. The Keldysh 
parameter has the value 1.87. 
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FIG. 2. Quantum mechanical (dotted line) and semiclassi- 
cal (solid line) ATI spectra for the Hamiltonian of Eq. ([!]) with 
Eo = 0.049a.u., loo = 0.09a.u. and the Gaussian potential Eq. 
&■ 

We obtain a quantum mechanical ATI spectrum (dot- 
ted line in Fig. g) with six distinct maxima. The semi- 
classical spectrum (solid line) is not shifted, the location 
of the maxima agrees with quantum mechanics. Hence, 
one can conclude that the softcore potential is respon- 
sible for the shift. The height of the third maximum is 
clearly underestimated and the details of the spectrum 
are exaggerated by the semiclassical calculation. Apart 
from these deviations the agreement is good enough to 
use this type of calculation as a basis for a semiclassical 
understanding of ATI. 

B. Semiclassical interpretation of the ATI spectrum 

1. Classification and coherence of trajectories 

With the chosen parameters most of the trajectories 
ionize during the pulse (~ 92%). We consider a trajec- 
tory as ionized if the energy of the atom 

s(t)=p(tf/2 + V(q(t)) (13) 

becomes positive at some time t n and remains positive, 
i.e. e(t) > for t > t n . Typically, the trajectories ion- 
ize around an extremum of the laser field. Tunnelling 
can not be very important, otherwise the agreement 
between quantum mechanics and semiclassics would be 
much worse. The Keldysh parameter of 1.87 suggests 
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that we are in between the tunnelling and the multipho- 
ton regime. Interestingly, the semiclassical description 
is successful although we are way below energies of the 
classically allowed over the barrier regime. 
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FIG. 3. Energy e(t) from Eq. ( |13| ) for trajectories ionized in 
the intervals Ji (solid line), I2 (dashed line), 73 (dashed-dotted 
line) and I4 (dotted line), respectively. For comparison, the 
laser field is plotted in arbitrary units (thick dashed line). 

An obvious criterion for the classification of the tra- 
jectories is the time interval of the laser cycle into which 
their individual ionization time t n falls, see Fig. ||. Typi- 
cally ionization of trajectory happens around t n = (2n — 
l)T/4 when the force induced by the laser reaches a max- 
imum. Hence, the ionized trajectories can be attached to 
time intervals /„ = \{n- l)T/2,nT/2], In Fig. |we have 
plotted four trajectories from the intervals I\ to Zj which 
end up with an energy E = 0.36 a. u.. After ionization 
each trajectory shows a quiver motion around a mean 
momentum pf Jl8| . One can distinguish two groups of in- 
tervals, namely those with trajectories ionized with pos- 
itive momentum pj (the intervals Izk-i) and those with 
trajectories with negative pj (the intervals fak)- These 
two groups contribute separately and incoherently to the 
energy spectrum as one might expect since the electrons 
are easily distinguishable. One can see this directly from 
the definition Eq. (^J) of the electron energy spectrum. 
For relative high energies fiuj the (short-range) potential 
may be neglected in the Hamiltonian Hq and we get 



OO 
OO 

Re Je lut {^f\e- ip2t/2h \^ f ) dt 


00 

j 5(u;-p 2 /2h) \^ f (p)\ 2 dp 

— OO 

( V f (-V2hui) 2 + ®f(V2hLj) ^ (h/2uj) 1/2 



Hence, to this approximation, the ATI spectrum is in- 
deed given by the incoherent sum of two terms belonging 
to different signs of the momenta of electrons ionized in 
different time intervals as described above. 

Figure Fig. ||(a) shows that Eq. (|l4|) is a good approx- 
imation. 
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FIG. 4. Upper panel (a): Semiclassical spectrum as an in- 
coherent sum ct+(l}) + a-(ui) (dashed-dotted line) compared 
with the full semiclassical spectrum (solid line). Lower panel 
(b): Semiclassical spectrum <j+{uj) , constructed with trajec- 
tories from the intervals I2, I4, Is and Is (dotted) compared to 
the incoherent sum 5+ of spectra that belong to the intervals 
I2 to 7g (solid line). 

Only for small to the spectra do not agree, where the 
kinetic energy is comparable with the (neglected) poten- 
tial energy. 

Quantum mechanically, all contributions from trajec- 
tories which lead to the same momentum pf of the elec- 
tron are indistinguishable and must be summed coher- 
ently. To double check that the interference from dif- 
ferent intervals /„ is responsible for the ATI peaks, we 
can artificially create a spectrum by an incoherent super- 
position cr+ = o"2 + 04 + ""6 + °~8 of contributions from 
trajectories ionized in the intervals I^j- This artificially 
incoherent sum (Fig. |](b)) shows similarity neither with 
cr+(w) nor with any kind of ATI spectrum. 



2. Classical signature of bound and continuum motion in the 
laser field 

The great advantage of an ab initio semiclassical de- 
scription lies in the possibility to make dynamical be- 
havior transparent based on classical trajectories, par- 
ticularly in the case of explicit time dependent problems 
where our intuition is not as well trained as in the case of 
conservative Hamiltonian systems. The classical quanti- 
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ties enter semiclassically mostly through the phase factor 
exp (i[S qp (t) - p{t)q{t))/K) = exp[i$/h} (15) 



which each trajectory contributes to the wave function 
Eq. (H|). Although the prefactor C qp (t) in Eq. (^) may be 
complex itself, the major contribution to the phase comes 
from the effective action $ in the exponent of Eq. (15). 
Figure |5| shows the energy s of the atom and the accu- 
mulated phase One can recognize a clear distinction 
between a quasi-free oscillation in the laser field after the 
ionization and the quasi-bound motion in the potential. 
The latter is characterized by an almost constant aver- 
aged bound energy (e(t)) (Fig. |(a)) of the individual 
trajectory giving rise to an averaged linear increase of the 
phase (Fig. ||(b)). After ionization the phase decreases 
linearly with an oscillatory modulation superimposed by 
the laser field. The almost linear increase of $ without 
strong modulation of the laser field during the bound mo- 
tion of the electron is remarkable, particularly looking at 
the laser induced modulations of the bound energy seen 
in Fig. ||(a). The averaged slope of the phase (positive 
for bound motion, negative for continuum motion) cor- 
responds via d^/dt = —E to an averaged energy. The 
behavior can be understood by a closer inspection of the 
action 



*(t) = S qp (t) - p{t)q{t) 

= f (2T-H-p(r)q(T)-q(T)p(T))dT-qp. (16) 
Jo 

Here, T — p 2 (t)/2 refers to the kinetic energy and H 
to the entire Hamiltonian of Eq. (Q) , the dot indicates a 
derivative with respect to time, and q = q(t = 0). With 
the help of Hamilton's equations and a little algebra $ 
from Eq. (16) can be simplified to 



$(«) 



e(r) - g(r) 



dV 
dq 



dr 



(17) 



where e is the atomic energy Eq. ©. With Eq. @ we 
can quantitatively explain the slope of $ in Fi g, p (b). For 
the low energies considered the potential Eq. ([12 ) can be 
approximated harmonically, 



V(q) « -Vq + V aq 2 



(18) 



Averaging $ over some time yields then w Vot, 

for any bound energy of a classical trajectory since for 
an oscillator averaged kinetic and potential energy are 
equal. Indeed, the numerical value for the positive slope 
in Fig. ||(b) is 0.6 a. u. in agreement with the value for Vq. 

For the ionized part of the trajectories we may assume 
that the potential vanishes. The corresponding solutions 
for electron momentum p(t) follows directly from Hamil- 
ton's equation p = —Eq sinwot, 



p(t) = — cos(uj t) + p, 

LOq 



(19) 



where p is the mean momentum. Without potential the 
phase from Eq. ([l7j ) reduces to $(t) = — Jp 2 (t)/2 dr 
and we obtain with Eq. ( |l9| ) 

*c(t) 

= ^-sin(2w t) - ^-sinujot- (U p +p 2 /2)t (20) 

2lo LU^ 

with the ponderomotive potential U p = E 2 /4ujq. We 
note in passing that Eq. (|2C|) is identical to the time de- 
pendent phase in the Volkov state (see the appendix). 
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FIG. 5. Part (a) shows the atomic energy e = p 2 /2 + V(q) 
as a function of time for three trajectories from the intervals 
I2 (dashed line), I4 (dotted line) and I§ (dashed-dotted line), 
part (b) shows the corresponding phases $(£). 



3. Semiclassical model for ATI 

The clear distinction between classical bound and con- 
tinuum motion in the laser field as demonstrated by Fig. |5| 
and illuminated in the last section, allows one to derive 
easily the peak positions of the ATI spectrum. Moreover, 
this distinction also supports the so called strong field ap- 
proximation (e.g. Jj^Q) where electron dynamics in the 
laser field is modelled by one bound state and the contin- 
uum. While this is postulated in |^] as an approximation 
and justified a posteriori by the results the correspond- 
ing approximation is suggested in the present context of a 
semiclassical analysis by the full classical dynamics, i.e., 
the behavior of the trajectories, as shown in |[ There, we 
have seen that each classical bound motion leads to the 
characteristic linear increase of the phase. If the entire 
phase space corresponding to the initial (ground state) 
wave function is probed with many trajectories of differ- 
ent energy, the dominant contribution will appear at the 
bound state energy which implies 



$&(*) » I P t, 



(21) 
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where I p is the ionization potential. The time for which 
a trajectory does not fall into one of the two classes, 
bound or continuum, is very short (Fig. |5|). Hence, we 
can approximately compose the true phase $ = + 
$ c . However, we don't know for an electron with mean 
momentum p when it was ionized. Hence, we have to 
sum over all trajectories with different ionization times 
r but equal final momentum p — pf which leads to the 
propagated wavefunction 

*/(*,p)~ / dTexp[i/h{<f> b (T)+$c{t)-<S> c (T))} 

where the phase A is given by 

A m „ = I p + U p + p 2 /2 - (n + 2m)hujQ . (23) 

From Eq. @ and Eq. @ follows that ATI peaks ap- 
pear at integer multiples nhuo of the laser frequency 
when 



— = nhuio - I p -U p . 



(24) 



One can also see from Eq. (|22|) that the ATI maxima be- 
come sharper with each optical cycle that supplies ioniz- 
ing trajectories. Of course, this effect is weakened by the 
spreading of the wavepacket hidden in the prefactor of 
each trajectory contribution (see Eq. (g)) not considered 
here. 

Trajectories that are ionized during different laser cy- 
cles accumulate a specific mean phase difference. The 
phase difference depends on the number k of laser cycles 
passed between the two ionization processes: 



V 



A*(p) = kT U p + *-^„ p 
The trajectories interfere constructively if 



A$(p) = 2ttZ 



1 , I 

2 P = k^ ~~ Ip ~ Up - 



(25) 



(26) 



If an energy spectrum is calculated exclusively with tra- 
jectories from two intervals separated by k cycles there 
should be additional maxima in the ATI spectrum with 
a distance fiuo/k. 

As a test for this semiclassical interpretation of the 
ATI mechanism we have calculated three spectra with 
trajectories where the mean time delay between ionizing 
events is given by At = T, At = 2T and At = 3T. For 
the spectrum Fig. ||(a) we have used exclusively trajec- 
tories from the intervals 1% and I4 (At = T). One can 
see broad maxima separated by Tiujq in energy. Trajec- 
tories from the intervals I2 and I§ (see Fig. f5|(b)) form 
a spectrum where the maxima are separated by fuvo/2 - 
as predicted for At — 2T. In analogy the separation for 



the ATI maxima in a spectrum with trajectories from the 
intervals I2 and Ig is given by ?iwo/3 (Abb. §(c)). The 
interference of trajectories ionized in many subsequent 
cycles suppresses the non-integer maxima according to 
Eq. (p3|). If the field strength is high enough the atom 
is completely ionized during the first cycle. The oppor- 
tunity for interference gets lost and we end up with an 
unstructured energy spectrum. 
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FIG. 6. Semiclassical spectra calculated with trajectories 
from the intervals I2 and I4 (a), I2 and (b), and h and Is 
(c). 

In an extreme semiclassical approximation we would 
have evaluated the integral in Eq. (^) by stationary 
phase. The condition 

d/dr[<P b ( T ) - $ c (r)] = I P + P 2 (t)/2 - (27) 

leads to complex ionization times t n whose real part is 
periodic and allows for two ionizing events per laser cycle, 
close to the extrema of the laser amplitude. The deriva- 
tion is simple but technical, therefore we don't carry it 
out explicitely here. However, it explains the observa- 
tion that ionization occurs close to the extrema of the 
laser field and it also makes contact with the tunnelling 
process often referred to in the literature since the com- 
plex time can be interpreted as tunnelling at a complex 
"transition" energy. 

Clearly our semiclassical analysis as described here 
supports the picture which has been sketched in pc[| 
interpreting a quantum calculation. The authors as- 
sume that wave packets are emitted every time the laser 
reaches an extremum. The interference of the different 
wave packets gives rise to the ATI peaks. 

In the following we will discuss the process of higher 
harmonic generation (HHG) which is closely related to 
ATI. In fact, the separation into a bound and continuum 
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part of the electron description is constitutive for HHG 
as well, the prominent features, such as cutoff and peak 
locations, can be derived from the same phase properties 
Eq. (||^) as for ATI. However, there is a characteristic 
difference, how these phases enter. 



IV. HIGH HARMONIC GENERATION 

First, we briefly recapitulate the findings of plf| , where 
we have calculated the harmonic spectrum with the soft- 
core potential Eq. ([ll] ) . With our choice of a = 2 the ion- 
ization potential is given by I p = 0.5 a. u.. The laser field 
has a strength E = 0.1 a.u., a frequency wo = 0.0378 
a.u. and a phase S = ir/2. The initial wave packet with a 
width of (3 = 0.05 a.u. is located at qp = E /loq = 70 a.u.. 
Note, that the cutoff energy Eq in such a symmetric laser 
scattering experiment is given by 



Ec — I P + 2U P . 



From the dipole acceleration (see Fig. 0) 

dV(x) 



d(t) = - 9{t) 



follows by Fourier transform 



dx 



*(*) 



<t(ll>) — / d(t) exp(iujt) dt 



(28) 



(29) 



(30) 



the harmonic power spectrum (see Fig. 
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FIG. 7. Quantum (a) and semiclassical (b) dipole acceler- 
ation of higher harmonics according to Eq. rt29|) . 
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FIG. 8. Quantum (a) and semiclassical (b) spectrum of 
higher harmonics according to Eq. (|3C|). 

Clearly, our semiclassical approach represents a good 
approximation. The dipole acceleration shows the char- 
acteristic feature that fast oscillations (which are respon- 
sible for the high harmonics in Fourier space) show only 
up after some time, here after t = T . This is the first 
time where trajectories are trapped. Trapping can only 
occur if (i) t n = nT/2, (ii) the trajectories reach a turning 
point (i.e. p(t n ) = 0), and (iii) if at this time the electron 
is close to the nucleus (q(t n ) « 0). The trapped trajecto- 
ries constitute a partially bound state which can interfere 
with the main part of the wavepacket (trajectories) still 
bouncing back and forward over the nucleus driven by 
the laser. The group of briefly bound (i.e. trapped or 
stranded) trajectories can be clearly identified, either by 
their small excursion in space (Fig. |9ja) or by the positive 
slope of their action (Fig. ^|b) as it was the case for ATI 
(compare with Fig. ^|) . By artificially discarding the ini- 
tial conditions in the semiclassical propagator which lead 
to trapped trajectories one can convincingly demonstrate 
that the plateau in HHG generation is a simple interfer- 
ence effect (2l]] . Here, we are interested firstly in linking 
ATI to HHG by using the same separation in bound and 
continuum parts of the dynamics already worked out for 
ATI. Secondly, we want to go one step further and con- 
struct a wavefunction based on this principle. 

Semiclassically, we have to look first at the phases of 
the observable. Therefore, we define a linear combination 
for the wavefunction from the respective phase factors for 
bound and continuum motion. Considering only terms 
in the exponent the harmonic spectrum Eq. (J3fj) reads 
simply 
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FIG. 9. Examples for direct (solid line), trapped (dotted 
line), and stranded (dashed line) trajectories, see text. 



<j(u>) ~ J dt exp(iwi) 

|exp (t# c (t)/ft) +cexv{i<$> b {t)/h)Y 



(31) 



where c ^ is a (so far) arbitrary constant. In princi- 
ple, c = c(t), however its change in time is much slower 
than that of the optical oscillations of the phases $(t), 
hence we may approximate c by a constant. The bound 
and continuum phases, ${, and $ c , are defined in Eq. 
( ^l| ) and Eq. (|2(]), respectively. For $ c we have p = 0, 
since this is the dominant contribution from the center 
of the wavepacket which was initially at rest. The re- 
sult is shown in Fig. [Io| Indeed, the plateau with the 
harmonics is generated, however, the initial exponential 
decrease is missing since we have neglected all prefac- 
tors of the semiclassical wavefunction which describe the 
dispersion of the wavepacket. 
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FIG. 10. Harmonic spectrum according to Eq. (pl[). 

Consequently, one can evaluate Eq. ( |3l| ) in stationary 
phase approximation. The integrand of Eq. (p|) becomes 



dt 



[Tuvt±($ b (t)-$ c (t))] = 



which happens at 

huj = 2U P sin 2 (ojt) + I p . 
From Eq. (133) we conclude the cut-off law 



2U V 



(32) 



(33) 



(34) 



as expected for laser assisted electron ion scattering . 
Using the same expansion into Bessel functions as in Eq. 
( ^2| ) we obtain for the spectrum Eq. 



dt exp 

oo 

E 



(Huj — Up — I p )t+ — — sin (2uot) 
2ujq 



dt u p —i p +2khwa)/}i j 



U, 



p 



2hui, 



(35) 



k— — oo ' 

Therefore, we see maxima in the harmonic spectrum for 



2kujQ . 



(36) 



We can go one step further and construct a full time- 
dependent wavefunction from this semiclassical approxi- 
mation, namely 

$(x,t) = 1^ c (x,t) + c^ {x) exp(itl p /h). (37) 

Here, ^a^x) exp(il p t/h) is the time dependent ground 
state wave function (without the laser field) and ^ s p(x, t) 
is a (semiclassical) wavepacket in the laser field but with- 
out potential. Calculating the dipole acceleration and 
the resulting harmonic spectrum with this wavefunction 
leads to a remarkably good approximation of the true 
quantum spectrum (compare Fig. ^ with Fig. |TT| ) . The 
dispersion of the wavepacket leads to the lower plateau 
compared to Fig. RG. 
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FIG. 11. Harmonic spectrum, generated from the wave- 
function Eq. @ with c = 0.025 and = 0.05 a.u.. 



V. CONCLUSIONS 
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A. Semiclassical comparison between ATI and HHG 

Clearly, the main structure such as the plateau, cutoff 
(HHG) and the occurrence of peaks and their separation 
in energy (ATI and HHG) is a property of the differ- 
ence of the classical time-dependent actions $f,(t) — § c (t) 
alone. 

However, the HHG power spectrum Eq. (^0|) is an inte- 
gral over all the time for which the electron wavepacket is 
exposed to the laser field. In contrast, the ATI spectrum 
is obtained in the long-time limit t — ► oo after the laser 
has been switched off. This difference may explain why 
the HHG results tend to be better than the ATI results 
semiclassically: Any semiclassical approximation (which 
is not exact) become worse for large times. 

A second point refers to the fact that the character- 
istic phase difference $&(£) — &c(t) appears already in 
the wavefunction Eq. ( p2| ) for ATI, while for HHG it oc- 
curs only in the expectation value Eq. (p9|). However, 
this difference is artificial, since the expectation value, 
or better its Fourier transform the power spectrum, is 
not the observable of higher harmonic radiation. The 
correct expression is the dipole-dipole correlation func- 
tion R which can be approximated as R oc |cr(oi;)| 2 under 
single atom conditions or in the case of an ensemble of 
independent atoms which radiate |l0| , |22| ] . Hence, in both 
cases, ATI and HHG, the peak structure appears already 
on the level of the quantum amplitude (or wavefunction) 
and is amplified in the true observable. 



B. Summary 

We have given a time-dependent fully semiclassical 
description of multiphoton processes. The prominent 
ATI and HHG features emerge naturally from proper- 
ties of the classical trajectories whose contributions to 
the semiclassical wavefunction interfere semiclassically. 
Any effect of this semiclassical interference can be double- 
checked by disregarding the phases. This leads (with the 
same trajectories) to a classical observable. As we have 
seen, to a good approximation the classical action for 
an individual trajectory can be composed of one part 
for the time the electron is bound (disregarding the laser 
field) and of another part $ c for the time the electron 
is in the continuum (disregarding the atomic potential). 
The relevant phase difference — $ c leads in both cases, 
ATI and HHG, to the prominent harmonic structures in 
terms of the laser energy Tiujq. Finally, we have been 
able to construct a simple wavefunction for higher har- 
monics generated in laser assisted scattering. Its key ele- 
ment is an explicitely time-dependent wavepacket of the 
electron under the influence of the laser field. Starting 
from an initial Gaussian distribution localized in space 
the wavepacket disperses in time providing the correct 
decrease of the intensity of the lower harmonics and in 
turn the correct height of the plateau. 
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APPENDIX: 

We want to calculate the semiclassical wave function 
of a free particle in a laser field according to Eq. (^). A 
particle in a laser field Vl(x, i) = Eq sin(wi) moves with 



E 

p(t) = p-\ cos(wt) = p+p(t) 



(Al) 



E 

q(t) = q+pt+ -^sm(ujt) = q + pt + q(t) (A2) 
The weight factor C qp (t) is given by 



(A3) 



For the phase factor S qp (t) — p{t)q(t) we get: 
S qp {t) - p{t)q{t) = - ^ sin(2wt) - U p t 

p 2 

~^-t-q{t)p-iP (A4) 

Evaluating Eq. (^) with the stationary phase approxima- 
tion, which is exact for quadratic potentials, leads to the 
condition that 



i ( p 2 7 f3 

f(q,p) =T xp(t) - —t-q(t)p qp q p 

n V I a a 



7 



7/3 



1 



2hf 



(A5) 



must have an extremum. With 



§£ = = 7 [ x -q(t)}-^ (q-qp)- ^-p (A6) 
oq ana 

g=0 = 7 [x-g(t)]*-^p 

+ i(x-pt-m- 1 9-^) (A7) 
S \ a a I 



we find 



q s 

Ps 



x — q{t) + ihfitqp 



1 + ihfit 
ih(3 



1 + ih/3t 



(A8) 

{x - q(t) - q ( ,) . (A9) 







After some algebra we arrive at the stationary exponent 



/(?.,„.) = jfoK*)- 2(1H f w) 



(x - <?(t) - g,a) 



(a; - - q^) 2 , 



2er(t) 

where a(t) is given by 

cr(t) = 1 + (3 2 h 2 t 2 



(A10) 



(All) 



The determinant of the second derivatives of / still has 
to be calculated. With 



d 2 f _ 7 4 + 2jf3 

dq 2 a 

d 2 f i 7 

dqdp h a 



we get 

( &f d 2 f \ 

dq 2 dq dp 

d 2 f d 2 f 



dp- 



2 = -t-* 



l 

h 2 a 



(A12) 



dct 



^-([l-i-fM/2] [l + i/3M]) 



\dpdq dp 2 J 



(A13) 

The factor 7 cancels as it should be and we are left with 

1/4 



exp 
exp 
exp 



1 



1 + ihflt 



p{t)x 

i h 2 /3 2 



sin(2wt) - U p t 



h 2a(t) 
t_ 

2a{t) 



(x - q(t) - qp) t 



(x - </(*) - qp)' 



(A14) 



This semiclassical time dependent wavepacket is quan- 
tum mechanically exact and corresponds to a superpo- 
sition of Volkov solutions according to a Gaussian dis- 
tribution at time t = 23 1. The fact that the semi- 
classical wavefunction is exact is a direct consequence 
of the Ehrenfest theorem which implies that interactions 
V ex x™, n = 0,1,2 have quantum mechanically exact 
semiclassical solutions. 
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